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The peak luminosity - peak energy correlation in GRBs 
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ABSTRACT 

We derive the peak luminosity-peak energy (iiso-.Epeak) correlation using 22 long Gamma- 
Ray Bursts (GRBs) with firm redshift measurements. We find that its slope is similar to the 
correlation between the time integrated isotropic emitted energy E lso and £ pea k (Amati et 
al. 2002). For the 15 GRBs in our sample with estimated jet opening angle we compute the 
collimation corrected peak luminosity L 7 , and find that it correlates with -E p0 ak- This has, 
however, a scatter larger than the correlation between -E pca k an d E~ (the time integrated emit- 
ted energy, corrected for collimation; Ghirlanda et al. 2004), which we ascribe to the fact that 
the opening angle is estimated through the global energetics. We have then selected a large 
sample of 442 GRBs with pseudo-redshifts, derived through the lag-luminosity relation, to 
test the existence of the Li so -£peak correlation. With this sample we also explore the possi- 
bility of a correlation between time resolved quantities, namely Lf and the peak energy at 

peak' 



the peak of emission 
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1 INTRODUCTION 

Several correlations have been identified among the intrinsic prop- 
erties of the (small) population of GRBs with measured redshifts 
z. In particular two spectral correlations have been recently dis- 
cussed in the literature: i) the "Amati correlation" between the en- 
ergy Speak where most of the emission is radiated and the total 
emitted energy (isotropic equivalent) Ei so (Amati et al. 2002 - A02 
hereafter; Lloyd & Ramirez-Ruiz 2002); ii) the "Ghirlanda corre- 
lation" between -E pca k and the collimation-corrected energy E^ 
(Ghirlanda, Ghisellini & Lazzati 2004 - GGL04 hereafter). It is im- 
portant to notice that these correlations refer to the time integrated 
spectral properties of GRBs. This is true for both Speak and for the 
spectral indices required to calculate the rest frame bolometric Ei so 
and Ey. 

However, time resolved spectral analysis of large samples of 
bursts (e.g. Ford et al. 1995; Preece et al. 2000; Ghirlanda, Celotti, 
Ghisellini 2002) have proved that the GRB spectrum evolves in 
time during the prompt emission phase. The spectral evolution is 
different among different GRBs (e.g. Ford et al. 1995) and not 
clearly linked to other GRB global parameters (e.g. duration, num- 
ber of peaks, peak flux). This spectral evolution may be revealing 
of the time variation of the parameters of the radiative process(es) 
acting in GRBs (e.g. Liang & Kargatis 1996) and/or of the rela- 
tivistic properties of the emitting outflow (e.g. Ryde & Petrosian 
2002). In order to understand the origin of such correlations it is 
thus compelling to determine whether they are representing global 
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energetics characteristics or they hold for and are dominated by the 
time resolved spectral properties, as expected if determined by the 
emission process(es). One obvious possibility is to test them against 
the peak luminosity, well defined for all bursts with known z (e.g. 
Liang, Dai & Wu 2004). 

This issue has been recently considered by Yonetoku et al. 
2004 (Y04, hereafter). With a sample of 12 GRBs of known z 
they found that -E pca k oc L°^. This correlation appeared to be 
tighter (but with similar slope) than the E pf!a ,y-Ei BO correlation, as 
originally found by A02. Note that the Y04 analysis adopts E pBa ^ 
and the spectral indices of the time integrated spectrum and not the 
spectral properties at the peak flux. 

In this Letter we first re-examine the -Bpcak-^iso correlation 
(i.e. the "Yonetoku correlation") with an enlarged sample of 22 
GRBs with spectroscopically measured z and published spectral 
properties. For 15 out of these 22 GRBs we have an estimate of 
their jet opening angle 9j (GGL04). We can thus calculate the col- 
limation corrected peak luminosity L- t and verify if there exists the 
equivalent of the Ghirlanda correlation - namely L 7 replacing £J 7 
(Sec. 2). Then we consider a much larger sample of 442 GRBs with 
z estimated through the lag-luminosity correlation (Band, Norris & 
Bonnel 2004 - BNB04 hereafter) to test if the Yonetoku correlation 
still holds for this whole sample (Sec. 3). In Sec. 4, by means of this 
same sample, we also study the relation between L? so and -Ep Cak , 
(i.e. using spectral parameters at the peak of the flux), to check 
whether this correlation is tighter than the Yonetoku one and we 
discuss the differences between the two. We find that the Ghirlanda 
correlation has a smaller scatter than the corresponding E pea k-Lj 
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Figure 1. Rest frame peak energy -E pC ak = ^p^ak z ) versus tne bolo- 
metric peak luminosity. Samples: GRBs with measured z listed in Tab. 1 
(blue symbols) - upper/lower limits are excluded except for the 2 X-ray 
Flashes (stars), shown for comparison with Fig. 1 of GGL04; 15 GRBs 
with a firm measure of the jet break time (from Tab. 2 of GGL04) and 
hence of the jet opening angle 9j (open yellow circles); same 15 GRBs 



once corrected for the (1 



j) collimation factor (red filled circles). 



Fits: best power-law fit to the S pC ak-^iso correlation (dashed blue line); 
best fit E, 



peak" 



J pcak 

-L 7 correlation (dashed red line); the Amati correlation from 
23 GRBs in GGL04 and GGF05 (long-dashed blue line) and the Ghirlanda 
correlation from GGL04 (dot-dashed blue line). 



Figure 2. Distribution of the scatter of data points around their best fit 
correlations. Filled histogram: 442 GRBs with pseudo z (black crosses in 
Fig. 3) around their best powerlaw fit (solid blue line in Fig. 3), and Gaus- 
sian fit to this distribution (solid black line). Red hatched histogram: scat- 
ter of the 22 GRBs (blue symbols is Fig. 1) around their best fit line (long 
dashed blue line in Fig. 1), and the best Gaussian fit (red solid line). Also re- 
ported are the distributions of the scatter of the 23 and 442 GRBs of GGF05 
with respect to their correlation in the -Epeak-^iso plane. 



correlation. We give an interpretation of this result in Sec. 5 and 
draw our conclusions in the final Sec. 6. 

In this paper we adopt a standard ACDM cosmology with 
= 0.7, O m = 0.3 and h = 0.7. 



2 THE Liso-Speak CORRELATION 

The bolometric 7-ray luminosity can be defined once the prompt 
emission spectrum and the redshift 2 of the source are known. 
GRB spectra are typically described by the Band function 
N(a, (3, Speak) (Band et al. 1993), parameterized by low and high 
energy power-laws (of photon indices a and f3, respectively) and 
by peak energy i5 pca k in the vF v representation. 

The burst emission varies on short timescales (e.g. Ramirez- 
Ruiz & Fenimore 1999) and no universal temporal profile describes 
the "zoology" of burst light curves (e.g. Norris et al. 1996). How- 
ever, in most cases, a dominating peak, with flux <J> integrated in 
the observed energy band, can be identified in the prompt emission 
light curve. The rest frame, bolometric (e.g. 1-10 4 keV), isotropic 
peak luminosity, including the redshift-energy band correction, fol- 
lows straightforwardly. In Tab. 1 we report the peak luminosities of 
the 29 GRBs examined by GGL04, computed assuming the time 
integrated spectrum of each GRB (as from Tab. 1 in GGL04). 
No published spectrum was found for GRB 011121 (detected by 
BeppoSAX) and this is no further considered. 

Speak versus Li S o for the these GRBs are shown in Fig. 1 
(blue symbols). We omit upper/lower limits except for the two X- 



ray Flashes with measured 2. Note that the underluminous GRB 
980425, associated with SN 1998bw, and GRB031203, associated 
with SN 20031w, are major outliers for both the Yonetoku and the 
Amati correlations. The statistical results for the correlations are re- 
ported in Tab. 2, together with the corresponding best fitting power- 
law parameters (weighting for the errors on both coordinates). The 
highly significant correlation has slope similar to that found by Y04 
with 12 GRBs. The distribution of the scatter measured along the 
correlation (i.e. the distances of the data points from the fitting line) 
of the 22 GRBs is shown in Fig. 2 (red-hatched histogram). A 
Gaussian fit (red solid line) to the distribution yields a scatter com- 
parable to that of the 23 GRB in the Speak-Siao plane (Ghirlanda, 
Ghisellini & Firmani (2005) - GGF05 - and black dashed line in 
Fig. 2). 

For 15 out of the 22 GRBs listed in Tab. 1 we can correct 
their isotropic luminosity Li so for the jet opening angle (9j (Tab. 2 
in GGL04), i.e. L 1 — Li so (l — cos 6j), with a corresponding error 
given by: 




The red symbols in Fig. 1 define the L 7 -_B poa k correlation. Again 
all the statistical parameters are reported in Tab. 2. The scatter of 
the best fit correlation (dashed red line in Fig. 1) decreases with 
respect to that using Li so - a trend similar to that found going from 
E iso to E 1 (GGL04 and GGF05). As discussed by GGF05 in rela- 
tion to the Amati correlation, also the scatter of the Yonetoku cor- 
relation found here can be interpreted as due to the distribution of 
jet opening angles. Note, however, that the scatter in the i5 pe ak-£ 7 
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Figure 3. Rest frame peak energy -E pca k versus peak luminosity Li so . 
Samples: 442 GRBs with pseudo z defined in GGF05 (black crosses); 22 
GRBs of Tab. 1 (red symbols). Fits: the powerlaw fit to the 442 black crosses 
(22 red symbols) is represented as solid blue line (dashed red line). The 
shaded regions show the lcr and 3<r width of the dispersion of the black 
crosses around their best fit. 



correlation is larger than that of the Ghirlanda correlation _E pca k- 



3 THE PSEUDO REDSHIFTS SAMPLE 

The original Amati correlation was found with 9 BeppoSAX GRBs 
with known z. Through a redshift independent test, Nakar & Piran 
(2004) and Band & Preece (2005) claimed that the larger BATSE 
sample is inconsistent (at 40% and 88% level, respectively) with 
the original Amati correlation. However, GGF05 (see also GGL04) 
have confirmed the above correlation (but finding a larger scatter) 
using a sample of 23 bursts with measured z as well as using a 
sample of hundreds GRBs with pseudo z. An even more general 
conclusion, i.e. the consistency of the above correlations with the 
entire BATSE long bursts sample, has been derived by Bosnjak et 
al. (2005). 

The same test of GGF05 can be performed for the -E pea k--Li S o 
correlation found in Sec. 2. More importantly, it is worth to inves- 
tigate if its scatter and slope change using this much larger sample. 
To this aim we consider the same sample defined in GGF05, which 
comprises 442 GRBs with pseudo z [estimated by BNB04 through 
the lag-luminosity relation] and known peak energy of the time 
integrated spectrum (found by Y04). Since the photon spectral in- 
dices are not given in Y04, in order to compute Li so we assume 
typical values, i.e. a — —0.8, f3 = —2.5 (see e.g. Preece et al. 
2000). 

In Fig. 3 we show these 442 GRBs (black crosses) in the rest 
frame -E pca k vs Li so plane, and in Fig. 2 we show the distribution of 
the scatter of the 442 points around this correlation (blue histogram 
in Fig. 2) together with its Gaussian fit (black solid line in Fig. 2). 
This scatter is only slightly larger than that defined by the 22 GRBs 



with measured z (red-hatched histogram and red solid line in Fig. 
2), and the slope of the correlations is similar (see Tab. 2). 

For comparison with what found for the Amati correlation by 
GGF05, Fig. 2 also reports the scatter distribution for the -E pea k- 
i?iso relation with the same samples of 23 and 442 GRBs (dot 
dashed red and black dashed line, respectively). As already men- 
tioned, the scatter of the Yonetoku correlation for the 442 GRBs 
can be interpreted as due to the distribution of the jet opening an- 
gles. Assuming that the _B pea k-L 7 correlation has a smaller scatter 
than the Yonetoku one, we can estimate the jet opening angle dis- 
tribution for the 442 GRBs: we find a lognormal with a peak at 
9 ~ 5°, i.e. consistent with that found in GGF05. 



4 THE Epeak-^Fso CORRELATION 

To the aim of investigating the spectral correlations at the peak of 
the prompt emission, the most correct approach would be to ana- 
lyze the spectrum of each GRB, time resolved at the burst peak. 
This would allow to derive a peak spectral energy -E pcak and a lu- 
minosity (from the spectrum at the peak of the burst) Lf Bo which, 
in general, might be different from the analogous integrated quan- 
tities. 

Given that the GRBs listed in Tab. 1 were detected by dif- 
ferent satellites and that the data are public only for BATSE, we 
can investigate the E^ ciik -Lf so correlation only with the sample of 
442 GRBs with pseudo z. For these, in fact, Mallozzi et al. (1998) 
provide the spectral parameters of the peak spectrum, derived by 
integrating the GRB signal for ~ 2 sec around the light curve peak, 
and BNB04 report the peak flux corresponding to the same peak 
spectrum. 

The sample of 442 GRBs allows a direct comparison with the 
results obtained in GGF05. However, we here exclude a few bursts 
because either their -B peak is below the BATSE ~ 30 keV energy 
threshold (7 cases) or -E pcak is not constrained by the spectral fit (1 1 
cases with /3 > —2). We report in Fig. 4 the remaining 424 GRBs 
with pseudo-jz. Also in this case we find a strong correlation with 
a scatter consistent with that obtained adopting Li so (Sec. 3), and a 
slightly flatter slope. In other words, the peak energy and luminos- 
ity at the peak are correlated, but the correlation is not significantly 
tighter than the Yonetoku correlation (see Tab. 2). 



5 THE ORIGIN OF THE SCATTER OF THE E PE ak-L-, 
CORRELATION 

In Sec. 2 we have derived the equivalent of the Ghirlanda correla- 
tion with the peak luminosity, i.e. E pea ^-L-, . This correlation (red 
symbols in Fig. 1) has a scatter a factor 1.7 larger than that of the 
Ghirlanda correlation (for the same GRBs). This implies that time 
integrated quantities correlate better than time resolved ("instanta- 
neous") ones at the peak of the emission. 

One possible reason can be envisaged by comparing the peak 
luminosity Li so as a function of the isotropic energy Ei so for the 
sample of 15 GRBs of measured z and jet opening angle (Fig. 5). 
The two quantities are correlated, but with a considerable scatter, 
which is comparable to the scatter of the _B poa k-L 7 correlation: in 
the insert of Fig. 5 the Gaussian fits to the scatter distributions of 
the Li so -_Ei so (a — 0.33, black line) and the -E pca k-L 7 correlation 
(a — 0.31, red line) are reported. Note that in this case the scat- 
ter corresponds to the 'horizontal' distance of the points from the 
fitting line. 




Figure 4. Rest frame peak energy E p , versus peak luminosity LP for 

° r aJ peak L J iso 

the sample of 424 GRBs with pseudo z. The red symbols represent the 22 
GRBs of Tab. 1 (reported for comparison although their peak energy and 
luminosity are derived assuming their time average spectrum, see Sec. 2). 
The solid blue line (dashed red line) is the powerlaw fit to the 442 black 
crosses (22 red symbols). The shaded regions represent the Icr and 3cr width 
of the dispersion of the black crosses around their best fit. 



Figure 5. Peak isotropic luminosity L lso versus isotropic energy i?i so for 
the 15 GRBs with measured z and jet angle. Different colors represent dif- 
ferent 2 ranges: z < 0.7 (blue), 0.7 < z < 1.55 (green), z > 1.55 (red). 
The solid line (with slope 1.0) shows the least square fit. The insert re- 
ports the Gaussian fit to the scatter of the data points around this correlation 
(black line), compared to the scatter (red line) of the 15 GRBs with known 
angle (red points in Fig. 1) around the Ep^^-Lj correlation. Also repre- 
sented (blue line) is the scatter of the same 15 points in the -Epcak-^iso 
plane around the Ghirlanda correlation. 



The smaller scatter of the Ghirlanda correlation with respect 
to the Epc^k-L-f one can be then ascribed to the fact that L~, has 
a larger spread with respect to E-y. Interestingly, this might be re- 
lated to E~j being a time integrated quantity, better representing the 
total kinetic energy of the fireball, which is the quantity involved in 
the derivation of the jet opening angle (e.g. Sari, Piran & Halpern 
1999). On the contrary, L 7 might be more subject to local and tem- 
poral fluctuations in the fireball, not representative of the global 
energetics. 



6 CONCLUSIONS 

We derived the -E pca k-£iso correlation with the current largest 
available sample of 22 GRBs with known z and well determined 
spectral properties (GGL04). This correlation has a slope 0.51, i.e. 
similar to that proposed by Y04, although its scatter is much larger 
than they originally found with 12 GRBs. The scatter is instead 
comparable with what GGF05 found for the Amati correlation us- 
ing the same sample of GRBs. Using the 442 GRBs with pseudo-z 
we still find a strong correlation, with similar scatter and slightly 
flatter slope than those found with the 22 GRB of measured z. 

We then considered the robustness of correlations for quanti- 
ties calculated at the peak of the emission, with respect to time in- 
tegrated properties. The former, in particular £7p Cak vs Lf Bo , results 
in a correlation equally tight to that involving integrated quantities. 

Correcting Li so for collimation, we found that the correspond- 
ing -Epcak-I^ correlation has a slope flatter than the Ghirlanda cor- 
relation (0.57 vs 0.7) and a larger scatter (0.17 vs 0.1). The larger 
scatter might be ascribed to the fact that the peak luminosity is less 
representative of the global energetics of the burst, which in turn is 



adopted to represent the total kinetic energy of the fireball and thus 
estimate the jet opening angle (e.g. Sari et al. 1999). 
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Table 1. Peak fluxes and bolometric luminosities for GRBs with measured 2 (in GGL04). S=BeppoSAX, B=Batse, I=Integral, H=Hete-II. Photon peak fluxes 
or energy peak fluxes with references (col. 2) and corresponding observed energy band (col. 5). References: (1) Amati et al. 2002; (2) Yonetoku et al. 2004; 
(3) 4fh Batse catalog; (4) Frontera et al. 2001; (5) Price et al. 2002; (6) Piro et al. 2004; (7) Sakamoto et al. 2004. (8) Sazonov et al. 2004. 
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Table 2. Statistical results. TV: Number of objects; r s : Spearman correlation coefficient; P: Chance Probability; r z : Partial correlation coefficient subtracting 
the effect of 2 (e.g. Wall & Jenkins 2003); A, S a and f2: Linear Fit Normalization, Scaling and Slope, i.e.: (E pe ak/100 keV) = A(S/ S ) S ; ft and a: Gaussian 
fit parameters. 
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